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Abstract. We apply a semiclassical approach to express finite temperature dynamical correlation functions 
of gapped spin models analytically. We show that the approach of [A. Rapp, G. Zarand, Phys. Rev. B 74, 
014433 (2006)] can also be used for the S = 1 antiferromagnetic Heisenberg chain, whose lineshape can 
be measured experimentally. We generalize our calculations to O(N) quantum spin models and the sine- 
Gordon model in one dimension, and show that in all these models, the finite temperature decay of certain 
correlation functions is characterized by the same universal semiclassical relaxation function. 



PACS. 75.10.Pq Spin chain models 
and statistics 



05.30.-d Quantum statistical mechanics - 05.50.+q Lattice theory 



1 Introduction 

The investigation of one-dimensional systems has a long 
history. In the past decades more and more one-dimensional 
and quasi-one-dimensional systems have been discovered 
and investigated experimentally. The discovery of these 
systems also triggered renewed interest in one-dimensional 
(1+1 dimensional) models and quantum field theories on 
• the theoretical side. In fact, many of the simplest though 
' non-trivial spin models and quantum field theories are re- 
alized in some of these experiments, which thus provides a 
unique possibility to compare theoretical predictions and 
reality. One-dimensional models provide the simplest ex- 
amples of spin liquid systems, can display spin-charge sep- 
aration, and they can also exhibit quantum phase transi- 
tions. Spin chains play a rather special role among one- 
dimensional systems. In 1931, Bethe solved the 5=1/2 
antiferromagnetic Heisenberg model, and showed that its 
spectrum is gapless [1]. His method has been general- 
ized, and is now widely used to study integrable one- 
dimensional systems. More than 20 years ago, however, 
Haldane pointed out that in contrast to the spin 5 = 
1/2 Heisenberg model, quantum Heisenberg chains with 
an integer spin do have a gap in the excitation spec- 
trum [2]. This astonishing difference between the half- 
integer and integer spin cases is by now accepted and 
well understood analytically [3] and numerically [4], and 
the existence of the gap was also confirmed experimen- 
tally [516171819110111] . Although zero temperature integer- 
spin chains and their static properties at finite tempera- 



tures have been studied extensively, understanding their 
finite temperature dynamical behavior still poses a ma- 
jor theoretical challenge. Not even the dynamics of the 
5 = 1 Heisenberg antifcrromagnct has been fully under- 
stood so far. Characterizing these dynamical properties in 
detail would be, however, crucial in order to interpret in- 
clastic neutron scattering [5|6|7|8|9] or nuclear magnetic 
resonance |10ll2j experiments. One of our goals in this pa- 
per shall be to describe the dynamical correlations of the 
5=1 antiferromagnetic Heisenberg chain, which are re- 
lated to the inelastic lineshape measured directly via neu- 
tron scattering by a Fourier transform. So far two methods 
proved to be efficient to study the dynamical correlations 
at finite temperature: 

(A) a field theoretical approach based upon a form 
factor expansion [13I14|15|16|17| of the dynamic structure 
factor (DSF). This method suggests that the lineshape of 
the dynamical susceptibility of the 5=1 antiferromag- 
netic Heisenberg chain around q = tt is rather asymmetric 
around the gap, and this asymmetry grows with increasing 
temperature |18j . 

(B) in an approach due to Sachdev and Young, the 
quasiparticles are described semiclassically and only their 
collisions are treated quantum mechanically [19] . 

In this paper we shall use the second, semiclassical 
method to calculate the finite temperature dynamical cor- 
relation function analytically for various spin models. This 
approach is based on the observation that in the spin mod- 
els investigated here, quasiparticles have a dispersion re- 
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Fig. 1. (Color online) The reflective S-matrix. During the 
semiclassical evolution of the system, the sequence of the in- 
ternal quantum numbers of the elementary excitations remains 
invariant. 



lation efc with a spectral gap A. Since the quasiparticlc 



density, p 



-A/T 



, is exponentially small at low tem- 



peratures, T -C A, the thermal de Broglie wavelength, 
AdB ~ 1/ VT, is negligible compared to the average separa- 
tion of the quasiparticles, ~ p~ 1 . Hence the quasiparticles 
behave quasiclassically except for the unavoidable colli- 
sions in d = 1 dimension. Generic gapped systems have 
quasiparticles which posses some internal quantum num- 
ber k. For a number of models [20I19I21I22I23] , it has been 
verified that at low temperatures, in the long wavelength 
limit the iS-matrix describing quasiparticlc collisions be- 
comes universally reflective (sec FigfTJ), 



(1) 



We shall use this simple structure of the (S-matrix to 
compute certain time-dependent correlation functions. In 
particular, we shall evaluate correlation functions of the 
form 

C(x,t) = (d(x,t)6(Q,0)) . (2) 

Here the operator 0(x, t) depends on the particular model 
studied: it corresponds to the operators S z (x,t) in the an- 
tiferromagnetic Heiscnbcrg chain and n z (x,t) in the 0(3) 
quantum rotor model, and it will be a vertex operator in 
the sine-Gordon model. 

Although the correlation function in equation ([2]) has 
already been computed for the sine-Gordon model by Damlc 
and Sachdev in reference |21| analytically, they determined 
the n z — n z correlation function for the 0(3) rotor model 
only numerically in their previous work, reference |22j . 
Here we show in detail that this correlation function can 
be obtained also analytically for all 0(N) models using 
the method of reference [5D] . As we shall see, for the mod- 
els considered here, the semiclassical relaxation function 
assumes the same universal form, as already stated in ref- 
erences [20] and [22]. 

In many cases, the operator O creates (or destroys) a 
quasiparticlc Therefore, at finite temperatures, we factor- 
ize the correlation function as 
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Fig. 2. (Color online) Relaxation function i? 7 (a;,t) for 7 = 
1/2. The inset shows the long time behavior i S> 1, the dashed 
lines correspond to the asymptotic diffusive form in equa- 
tion (JU). 



where C(x, t)r=o corresponds to the coherent propagation 
of a quasiparticlc, while the relaxation function i? 7 (x,i), 
defined through equation ([3]), describes the scattering with 
the thermally excited quasiparticles. For all models con- 
sidered here in the semiclassical limit, the relaxation func- 
tion is given by 



Ry(x, t) 



d(j) (1 — 7 2 ) cos (sin(0)x) 
2-7T 7 2 + 27 cos 0+1 



(4) 



exp 



1*1(1 — cos 0) \-^e- u + uerf(u) 

/7T 



where 7 is a model-specific parameter, and u = x/i, x = 
x/£ c and t = t/r c are the dimensionless velocity, position 
and time. Here the characteristic length and time scales 

z A / T ' and t, 



are defined as £ c = 7 
respectively, with A = 



2tt c 2 



T A ' c 1 1 

e(fc) I fc >o the quasiparticle gap and 

c 2 = A[d 2 ek / dk 2 ]\k^o the square of the characteristic "ve- 
locity of light" . The relaxation function is shown in Fig. [2] 
and Fig. [3j for the parameters 7 = 1/2 and 7 = 1/3, 
respectively. In the Q-state quantum Potts model, 7 = 
1/(0-1) EQ]. We shall prove that 7 = 1/N for the 
O(N) quantum rotor models, while in the sine-Gordon 

model we find 7 = —cos O^p-j > with 7' and 77 defined 
later [2T]. Note that in the case of the sine-Gordon model, 

1 \£j_ e A/T ^ respectively. 
The relaxation function for 7 < 1 is diffusive for long time 
scales 1<£ 



£c — \ y^TA 1 and t c - 



FLyix, t) 



1 



(5) 



where D 



C(x,t) = C{x,t) T =oRy(x,t), 



(3) 



tti= is the dimensionless diffusion constant. 

2y/7T 

On the other hand, in the transverse field Ising model, 
7=1 (Q = 2). In this special limit, 7 — > 1, the relaxation 
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Fig. 3. (Color online) Relaxation function _R 7 (a;,t) for 7 = 
1/3. The inset shows the long time behavior i ^> 1, the dashed 
lines correspond to the asymptotic diffusive form in equa- 
tion (|5|. 



becomes exponential instead of being diffusive |19j . Ex- 
ponential behavior also appears under the assumption of 

a transmissive/non- interacting scattering matrix — ¥ 

(— 1)5kISk2 HO]- We believe, however, that generic lattice 
spin models always have a reflective scattering matrix, and 
special conditions and cut-off schemes are required to ob- 
tain a diagonal 5-matrix. 

This paper is organized as follows. We shall discuss 
the semiclassical correlations of S = 1 antifcrromagnetic 
Hcisenbcrg model and their relation to the correlation 
functions in the 0(3) quantum rotor model in Sec. [5J 
These results will be generalized to the O(N) rotor chains 
in Sec. [3] In Scc.|4] we derive using our approach the semi- 
classical relaxation function of the sine-Gordon model, 
first derived in a somewhat different way and form by 
Damlc and Sachdev [21] . The conclusions are presented in 
Sec. [5] 



2 The AF Heisenberg chain and the quantum 
rotor model 

The 5=1 Heisenberg chain is defined by the Hamiltonian 
tf Hcis = J^SiSi+i , (6) 



where Sj is an S = 1 spin at site i, and J > is the 
antiferromagnetic exchange coupling. As first proved by 
Haldane, antiferromagnetic spin waves have a gap, if 5* 
is an integer [2j. Haldane's prediction is in good agree- 
ment with many experiments on different quasi-lD ma- 
terials |5I6I7I8I11I10I9] , although the SU(2) symmetrical 
Heisenberg model provides usually only an approximate 
description of these experimental systems, since pertur- 
bations like magnetic anisotropy, lattice distortions, dis- 
order, etc., are always present in physical systems. 



Low-energy excitations of the integer spin Heisenberg 
model in equation ([6]) can be mapped onto the 0(3) rotor 
chain }2I3| . defined by the following Hamiltonian, 



H 



J 9 



(T) 



Here hi is the position operator (orientation) of the rotor 
on site i with the constraint nf = 1, and Li = x f>i is its 
angular momentum operator satisfying the commutation 
relations [L^,L^] = ie a pryLl and [Lf,n^] = iea^nj. The 
mapping can be done most straightforwardly in the path 
integral formalism by representing spins in terms of coher- 
ent states N;(t) Si(r). The unit modulus field N;(t) 
can then be parametrized [3|2|24j by slowly varying fields 
as 



Ni(r) = (-l) 4 n,(r)^l - a 2 L?(r) + aL,(r) , (8) 

where n describes the staggered and L the uniform compo- 
nent of the Heisenberg spins, and a is the lattice spacing. 
In the S = 1 case, the excitations corresponding to both 
fields n(r) and L(r) are massive [2]. Integration over L 
leads to the non-linear u-model, defined by the action 



P 

A nlam = -J-- J dr J dx[(d T h( 



x,t)) 2 +5 2 (d x h{x,T))< 



(9) 

where g is a coupling constant, c has a dimension of veloc- 
ity, and the field h(x, t) satisfies the constraint n 2 (x, r) = 
1. The non-linear a- model is the continuum theory of the 
rotor chain, equation ([7]) [24] . 

Let us now discuss the properties of the rotor model 
defined by equation j7|. For g >• 1, in the ground state, 
all rotors need to be in the L = state to minimize kinetic 
energy. The lowest energy excitations form a triplet with 
quantum numbers L z = A = — 1, 0, 1. In d = 1 dimension, 
the qualitative structure of the low-energy spectrum is 
the same for any g > and the excitations have a gap 
A(g) |22|24] , At low enough temperatures, T < A, the 
gap allows us to apply a semiclassical approximation. 

Let us focus on the dynamical correlation function of 
the rotor model, and follow the same steps as for the quan- 
tum Potts model [20 . In the semiclassical limit, the cor- 
relation function can be approximated as: 



C mtm (x,t) = (n z (x,t)h* (0,0))^ rotor 

{A„} V V 



(10) 



x {{x v , V v: \ u }\h z (x, t)h z (0, 0)\{x v , v v: \ v }) 

where the function P({x v , v v , A„}) is the probability den- 
sity of having quasiparticles with velocities v v and internal 
quantum numbers L z v = A„(A„ = ±1,0) at positions x v . 
In the low density limit, P({x„, A;,}) factorizes as 



P{{x v ,v v ,X v }) = 



1 1 



(11) 
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Fig. 4. (Color online) Semiclassical dynamics in the 0(3) 
rotor model. Due to the reflective structure of the S-matrix, 
the sequence of the quantum numbers A„ does not change in 
time. 

with the distribution of the velocities given by the Maxwell- 
Boltzmann statistics, 

Calculating the matrix elements for n z cos ff) with the 
first few spherical harmonics, one finds that n z either cre- 
ates a quasiparticle (L = 1) with L z = \ = at x = 
with some velocity v or destroys one already present in 
the configuration {x u , v v , The probability of latter is 
exponentially small since at T <gC A the quasiparticle den- 
sity p is low, so we shall neglect this. Due to the collisions 
with the thermally excited particles, there are only cer- 
tain configurations where the quantummechanical overlap 
in equation ([TO]) will be non-zero. Similar to the Potts 
model discussed in reference [3D], the 0(3) rotor model 
has a purely reflective scattering matrix |22l23j , 

S%& = i-^xl < ■ (13) 

This structure implies that the sequence of the quantum 
numbers A„ does not change in time (see Fig. [3J . Thus if 
the extra particle created by n z collides with N + parti- 
cles from the right and particles from the left, then 
to have a non-vanishing overlap, the first n = N + — iV_ 
particles to the right must also have A„ = 0. At time t a 
quasiparticle must be destroyed by n z (x, t) such that the 
overlap in equation (fT0|) docs not vanish. These observa- 
tions are analogous to the arguments given for the para- 
magnetic side of the quantum Potts model [2D]. Taking 
also the phase factors into account, we get an expression 
similar to that obtained for the paramagnetic side of the 
quantum Potts model 

C rotOT (x, t) = Gjffi{x, t)R{x, t) , (f 4) 

where the first term is the average of a phase factor coming 
from the propagation of a single particle from (0, 0) to 



(x,t) [22] 

C r T t°j(x, t) = |-AV^7c 2 - t 2 ) , (15) 

where Ko is the modified Bessel function of the second 
kind and Z is a non-universal quasiparticle residue. Well 
within the light cone, x <C ct, one can recover the Feynman 
propagator of a particle with a mass of A: 

^ 1<d|f '~ e " , "&(^)' (16) 
The second term in equation fTJJ is given by 

°° (— 1)" / \ 

R(x,t) = 22 -^H~ \ S n,^2j0(x-x v -v v t)-0(-x v )]/ 

71 — — OO U ' uS 

= RlMx/Sc,t/T e ). (17) 

In the middle expression, we already took the average over 
the angular momentum components {A„}, leading to the 
factor 1/3K The expression in the average is the proba- 
bility that the excited particle drifts to the right by n. The 
factor of (—1)" corresponds to the net phase accumulated 
in course of the scattering with the quasiparticles. The 
average over the velocities and coordinates can be carried 
out analytically to obtain the result on the right hand side 
of equation §F7§ [2D] . 

Having derived the semiclassical correlation function 
of the 0(3) rotor model, the correlation function of the 
antifcrromagnctic Hcisenberg chain 

C(x,t) = 0'(x,t)S'(O,O)), (18) 

can be approximated at low temperatures as follows. As it 
can be seen from equation JS} , and proved more rigorously 
in reference [21], while the fluctuations around q = in 
the Heisenberg chain are described by the L z — L z corre- 
lation function of the 0(3) quantum rotor chain, [23] the 
dynamics at q = it / a is related to the n z — n z correlation 
function. As a consequence, in the low temperature limit, 
the dynamic structure factors of the Heisenberg and rotor 
models are simply related as 

S zz (q = n/a + k,Lo)^ 5 rotor (fc, w) . (19) 

The correlation function that we calculated for the ro- 
tor model analytically, has been obtained by Damlc and 
Sachdev numerically using a stochastic sampling. The L z — 
L z correlation function has been evaluated analytically in 
reference 1221 . 



3 The O(N) rotor chain 

As a generalization of the 0(3) quantum rotor chain, let 
us now investigate the correlations of rotor models with 
O(N) symmetry, defined by the Hamiltonian 

^W = ^£L?-j£nA +1 , (20) 

i i 
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where is the N > 3-dimensional position operator of 
the rotor at site i with the constraint nf = 1. The angu- 
lar momentum operator components are given by L a p = 
n a pp — h^p a , where [n a ,pp] = i8 a p, and the square of the 
angular momentum operator is defined as 



(21) 



It is easy to see that the operator n a creates a member 
of an A^-componcnt quasiparticlc multiplct. The reason- 
ing is based on the Wigner-Eckart theorem: The ground 
state is structureless, and transforms as an 0(N) singlet 
under O(N) rotations. However, from the commutation 
relations it trivially follows that n a -s form the compo- 
nents of an .ZV-dimensional irreducible tensor operator. 
As a consequence, to have a non-vanishing matrix ele- 
ment, (/3\h a \0), the state (5 must transform according to 
the same irreducible representation as h a , moreover, the 
matrix elements must be proportional to 5 a ,0- This im- 
plies that single particle excitations of momentum k and 
having a finite overlap with n"|0) must form Af-fold de- 
generate multiplcts, \k,a)(a = 1,..,N). Moreover, since 
they transform as the components of n a , they must also 
be orthogonal to each other, and have the matrix ele- 
ments (k, P\h a \0) = A{k)5 a .f3, with A(k) a momentum- 
dependent form factor. 

Similar to the 0(3) rotor model, the continuum limit 
of the 0(N) rotor chain maps onto the 1+1-dimcnsional 
O(N) sigma model. The asymptotic 5-matrix of the con- 
tinuum O(N) non-linear a-model is purely reflective, and 
here we shall assume that this holds also for the lattice- 
regularized Hamiltonian |23j . 

It is thus straightforward to generalize our discussion 
of the 0(3) rotor model to this case and we find that, 
in the semiclassical limit, the correlation functions of the 
O(N) rotor chains can also be factorized as 



O(N) 



C 



where C^^\x, t) corresponds to the T = temperature 
propagation of a quasiparticle. Hence the universal relax- 
ation function R^(x,t) describes the decay of correlations 
of O(N) rotor and non-linear sigma models in the semi- 
classical limit. 



(x,t) = (h N (x,t)n N (0,0)) 



C%W(x,t)Ry= 1/N (x,i), (22) 



4 The sine-Gordon model 

As a next example, let us discuss the sine-Gordon model, 
defined as 



16tt 



l/T 

dr I dx 



(a c <Z>) 2 + -(d T <?) 2 - 5 2 cos( 7 '<2>) 



(23) 





m ? — + 


m 4 = - m 6 = - 










My = - 






m 5 = + 




nil = + 








P>1 
















X, 


x 2 x 3 


x 4 x 5 x 6 


x 7 



Fig. 5. A configuration of kinks and antikinks in the sine- 
Gordon model at t = 0. 



Here <P is the sine-Gordon field, g and 7' are coupling 
constants, and c has the dimension of velocity. Diffusive 
behavior of the sine-Gordon model in the semiclassical 
limit was found first by Damle and Sachdcv [21] . In this 
model the correlation function is defined by 



C$(x, t) 



( e «j*(*,t) e -iu*(o,o)^ iG 



(24) 



Now we shall prove that the relaxation of correlations in 
equation (j2"4")l is described by equation ([!]), with the pa- 
rameter 7 set to 



/27T7A 

T = ~ C0S {— ) 



(25) 



Within the semiclassical approximation and using the 
results of the previous subsections, it is easy to show the 
equivalency of the correlators of the sine-Gordon model 
and the Potts model. Let us follow here a slightly differ- 
ent path than the one in reference [21]. It is known that 
low-energy configurations of the sine-Gordon field can be 
described in terms of domains, separated by soliton and 
anti-soliton trajectories. These domains can be labeled by 
integers so that the field is approximately 



${x, t)~^j^ m v O{x - x v (t)) 



7 



(26) 



where x v (t) are the trajectories of the solitons and m v = 
±1 are the charges of the solitons (see Fig. [5]). The collision 
of the solitons is described by a reflective 5-matrix [21] . 

Within the semiclassical approximation, we find that 
the average over the kinetic variables and soliton charges 
factorizes as 



C${x,t)- £ {S n j2je( x - x „- Vu t)-e(o-x„))/ {xi/ ^ } 
x(-l)» /^P i (n)e ,,, 7( a -W)\ . (27) 

\ k I {m„} 

Here Pk{n) is the probability that the total number of 
m v — +l-charge soliton lines crossing the section (0,0) — ► 
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(x,t) is k, while the total number of collisions along the 
section is n (25]. The second term in equation ([27]) can 
be simplified as follows. The probabilities that a soliton 
or an anti-soliton is associated with a given trajectory 
are equal in a configuration. Therefore the probability for 
intersecting k soliton lines has a binomial distribution 

since every sequence of a given number of solitons and 
anti-solitons has equal probability. Using equation (|2"8"|) . 
we can compute the sum over k in equation (|27[) to obtain 
the result 

/^P fc (n)e^( 2fe -H)\ =cos L^L\ lnl - (29) 
\ fe / {m„} V 7 / 

Comparing this expression with the relaxation function 
found for the Potts model one finds equation (|2"5j) . 



5 Conclusions 

In this paper we obtained an analytical expression for the 
relaxation of certain correlation functions in a number of 
gapped one-dimensional systems, using the same scmiclas- 
sical approach we applied earlier for the Q-state quantum 
Potts model [5D]. In particular, here we studied the O(N) 
rotor models, discussed antiferromagnetic correlations in 
the S = 1 spin Heisenberg model, and also analyzed the 
dynamical correlations of the sine-Gordon model. 

Similar to the Potts model, all these systems have a 
gap in their quasiparticle spectra, and their quasiparti- 
cles posses some internal quantum numbers. In all these 
models, the collision of the quasiparticles is described by 
a reflective 5-matrix. The main result of our work is that 
the relaxation at finite temperature is given by the same 
universal function R^(x,t), defined by equation More- 
over, as shown in reference [20], the semiclassical relax- 
ation function i? 7 (x,t) also describes relaxations in the 
transverse field Ising model, corresponding to 7 = 1. 

For models with the parameter 7 < 1 ( spin models 
having a multiplct of excitations), the expression derived 
here predicts that the long time behavior is diffusive. In 
the case of the Q = 3 quantum Potts model on the ferro- 
magnetic side, the origin of this diffusive behavior is the 
random diffusion of the domain containing the origin at 
time t — 0, due to the random motion of the kinks which 
separate neighboring domains. On the other hand, in the 
case 7 — ► 1, the singularity at <fr = tt dominates the in- 
tegral in equation leading to an exponential rather 
than diffusive decay. This exponential behavior is a conse- 
quence of a "destructive interference ": for the Ising model 
in a transverse field, the sequence of domains on the fer- 
romagnetic side is always alternating. Carrying out the 
calculations assuming a non-interacting/transmissive 5- 
matrix also leads to exponential decay in the correlation 
functions [20] . 



We thus find that the relaxation function equation ([4]) 
is rather universal and it appears in essentially all one - 
dimensional models having a reflective 5-matrix. The idea 
of universal diffusive or exponential relaxations depending 
on the structure of the scattering matrix appeared already 
in references [20121] , and in reference [17] , although the re- 
sults of the latter work do not entirely agree with the relax- 
ation function obtained within the semiclassical approach. 
Apparently, the form-factor expansion of reference [17] 
builds in long-ranged correlations to the correlation func- 
tion, since at t = 0, it decays as ~ 1 5r | 1 / 2 . According 
to physical intuition, any spatial correlation between kink 
quantum numbers should decay beyond £ c . Therefore, on 
simple physical grounds, one expects exponential decrease 
in space for the equal time correlation function, which con- 
tradicts the result of reference Q2]. The It]' 1 / 2 behavior 
for x = 0, on the other hand, can be understood simply as 
a result of diffusive motion of the quasiparticles. [20] The 
authors of reference [T7] have also shown that for a diago- 
nal iS-matrix the decay becomes exponential. We believe, 
however, that the generic 5-matrix is of an exchange form. 
In fact, Damlc and Sachdev argued in reference [3T] that a 
transmissive 5-matrix needs the fine tuning of an infinitely 
large number of coupling constants. A similar observation 
was made in the discussion of the 5-matrix of the Q — 3 
state Potts model, where we argued that on a lattice, the 
scattering matrix is reflective in the low velocity limit. 

We are thus tempted to call the relaxation function 
Ry(x, t) the universal semiclassical relaxation function for 
gapped spin systems. We must mention that the diffu- 
sive behavior discussed here derives from the asymptotic 
form of the 5-matrix, and it therefore has a somewhat 
limited range of validity: At finite temperatures, the mo- 
menta of incoming particles remain finite, and therefore 
the scattering matrix is not purely reflective. We esti- 
mate that the probability that a quantum number flips 
rather than being reflected during a collision is Pfn p ~ 
a 2 {Ak 2 ) ~ a 2 TA/c 2 , where a is the lattice spacing and 
Ak is the momentum difference between the quasipar- 
ticles. These estimates carry over to all models consid- 
ered here, since in all these models the amplitude of non- 
rcflcctivc scattering in the leading order is proportional to 
the typical quasiparticle momentum. [20I21I22I23I26] As a 
consequence, the diffusive decay must be asymptotically 
replaced by an exponential decay at very long time scales 
and very large separations. The timescale at which the 
diffusive behavior is expected to vanish is tdiff ~ T cPf^p ~ 
(c 2 A/a 2 )e A / T . The corresponding lengthscale is Xdis ~ 
VDt^~c 2 /(aTA)e A / T . 
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